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ABSTRACT 

We calculate the thermal and dynamical evolution of the surface layers of an accreting neutron star 
during the rise of a superburst. For the first few hours following unstable ^^C ignition, the nuclear 
energy release is transported by convection. However, as the base temperature rises, the heating time 
becomes shorter than the eddy turnover time and convection becomes inefficient. This results in a 
hydrodynamic nuclear runaway, in which the heating time becomes shorter than the local dynamical 
time. Such hydrodynamic burning can drive shock waves into the surrounding layers and may be the 
trigger for the normal X-ray burst found to immediately precede the onset of the superburst in both 
cases where the Rossi X-Ray Timing Explorer was observing. 

Subject headings: accretion, accretion disks — nuclear reactions, nucleosynthesis, abundances — stars: 
neutron — X-rays: bursts 



1. INTRODUCTION 

Superbursts are hour-long X-ray bursts powered by un- 
stable thermonuclear burning on the surfaces of accret- 
ing neutron stars in low mass X-ray binaries. Discovered 
with long term monitoring campaigns by BeppoSAX and 
the Rossi X-Ray Timing Explorer (RXTE), their recur- 
rence time (« 1 yr), duration (« hours), and energies 
(« 10**^ ergs), are ~ 1000 tim es greater than that of nor- 
mal t ype I X-ray bursts (se e Kuulk er^ 120041: iCumminel 
120041: IStrohmaver &: Bildsten .2006 for reviews) . 

The current view is that superbursts are fu- 
eled by the u nstable burning of a thick layer of 
J^^C llCumming & Bildsten 200^ IStrohmaver fc Brownl 
l2002fl . which calculations show i s present in the ashes of 
H/He burning in normal bursts l [Bmgfii fc Bildstenlll998t 
iSchatz et al ] l2008at iCooner et al.ll200fi|) . Fits to suner- 
burst lightcurves suggest ignition column depths of 0.5 — 
3 X 10^^ g cm~^ and energy releases of 2 x 10^^ ergs g" 



implying ^^C mass fractions X12 > 10% (jCumming et al I 

I200fij) . For accretion rates M < 0.3MEdd (where MEdd is 
the global Eddington accretion rate), a larger ^^C frac- 
tion is needed {X12 > 20%) if conditions for the thermal 
instability are reached befor e the ^^ C bu rns stabl y away 
(jCumming fc Bildstedl200lt ICumming etalMOO^ . 

Previous theoretical studies have focused on two 
aspects of the ^^C burning scenario: the thermal 
structure of the surface layers just before ignition 
JCumming & Bildsten 2001; Strohmavcr & Brown 2002; 
[grown 2004; Cooper fc NaravanI 120051: ICumming et all 
I2OO6; .Cooper et al...200fi|) and the thermal evolution of 
the fully burned layers as they cool following the burst 
ICumming & Bildsten"'2001'; 'Strohmavcr fc Brown"2002| 
Eymming & Macbeth 2004; Cooper & Naravan 2005; 
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ICumming et al .''2006'). In this paper, we investigate the 
rise of the superburst starting from ^^C ignition. Unlike 
the cooling studies, which assume the entire layer of fuel 
burns instantly, we resolve the thermal and dynamical 
evolution during the burst rise. Understanding the de- 
tails of the rise is important because it sets the initial 
conditions for the cooling models. 

The principle result of our study is that burning is so 
explosive during a superburst rise that the heating time 
th = {dhiTb/dt)^^ becomes shorter than the dynamical 
time td = h/cs ~ 10~^ s, where Ti,, h, and Cg are the 
temperature, pressure scale height, and sound speed at 
the base of the burning layer. Such a small th is possible 
because the burning takes place in a strongly degener- 
ate la yer (cf. ^ooslcv fc Taam 1976; Taam & Picklui3 
Il978j) . During a normal burst, H/He ignites in a weakly 
degenerate layer and the burning terminates once radia- 
tion pressure dominates a.t Th ^ 2 x 10^ K; thro ughout 
the burn th > lO'^ s (e.g.. iWeinberg et alJl2006j) . 

Once th ^ td, convection cannot transport the nu- 
clear energy release and the base undergoes a local ther- 
monuclear runaway, rapidly burning all the ^^C within a 
thin shell. A hydrodynamic combustion wave forms and 
propagates into the surrounding fuel, likely driving shock 
waves into the overlying layers. By contrast, in normal 
X-ray bursts the energy is transported exclusively by con- 
vection and heat diffusion (see Bildsten 1998). 

Our primary aim here is to describe those ignition con- 
ditions that result in a hydrodynamic nuclear runaway. 
That is, we examine the range of parameters for which 
the minimum heating time satisfies th,min < td- We 
therefore limit our present analysis to the sequence of 
events leading up to the formation of a combustion wave. 
In §12 we describe the initial conditions of our rise models 
and in §|3we present our main results. We conclude in 
§ 01 and briefly discuss how the hydrodynamic nature of 
the burning can influence the observed rise. 

2. INITIAL CONDITIONS 

We investigate how the superburst rise depends on 
three ignition parameters: the column depth yi, at the 
base of the ^^C layer, the initial ^^C mass fraction 
X12, and the composition of the heavy ashes of H/He 
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burning. We assume a neutron star mass and radius 
oi M — IAMq and R — 10 km and adopt a plane- 
parallel approximation with a constant surface gravity 
g = (GM/i?2)(l + z) = 2.4 X 10" cm s'^ and redshift 
1 + z = 1.3. We assume electrons, ions, and photons 
supply the pressure and calculate the equation of state, 
volumetric ne utrino emissivi ty e^, and thermal conduc- 
tivity K, as in lBrownl l)2004j) . 

The thermal profile of the accumulated layer is found 
by integrating the entropy and heat equations 



dy 

dT _ F 
dy pK ' 



(1) 
(2) 



where y is the rest mass column depth, F the heat flux, T 
the temperature, p the density, and Cnuc the rate of heat- 
ing from nuclear reactions. When the accreted layers are 
in hydrostatic balance, the pressure/? = gy. We integrate 
equations 1^ and (0) inwards, setting the outer boundary 
aX y = 10^ g cm~^, where we adopt a radiative zero so- 
lution. The flux at the outer boundary has contributions 
from hot CNO burning i^cNO, an d deep crustal heatin; 
Fnni.st, - We calculate Fcno as in iCumming fc Bildste 
lj2000(l and assume the accreted gas has a solar composi- 
tion. 

For a given choice of y;,, X12, and ash composition, 
we iterate to find the value of i^crust that results in un- 
stable -'^^C ignition at y^. We define unstable ^^C igni- 
tion according to a (numerically solved) one-zone thermal 
instability criteri on dt„ur.ldh\T = dtr. nr,\ / d\n.T , where 
e^nni = oKTlv'^ (iFuiimoto et alJ I lQSlt iFushiki fc LambI 



Il987t ICumming fc Bildstenll20nni: [Coming et alJl20na) 

We change the composition from H/He fuel to ashes 
of H/He burning at the depth j/hc where He unsta- 
bly ignites^ assuming a local accretion rate per unit 
area m = O.lmEdd, where WEdd is the local Edding- 
ton accr etion rate. The ^^C energ y generation rate is 
given bylCaugh lan fc F owled l|1988(l with screening from 
IQgata et al.l 1)19931. This rate only accounts for the en- 
ergy released in fusing ^^C to ^^Mg and is thus an un- 
derestimate if the ^^C burns to iron-group elements or if 
the heavy (trans-iron) ashes of H/He burning p hotodis- 
integrate ijSchatz. Bildsten. fc Cummindl2003b(l . 

3. THE SUPERBURST RISE 

The rise of a superburst can be divided into three 
stages: a convective stage, a runaway stage, and a hydro- 
dynamic stage. During the convective stage, convection 
transports the heat flux from the base yh to the top of the 
convective zone j/c — Vc{Tb)- We calculate ydTb) in S I3.1I 
and show that since the thermal time tth(2/c) = CpVil pK 
is greater than the heating time th CpTb/ e^uc, the con- 
vective zone extends outward to lower pressures as Tf, 
rises (Cp is the specific heat at constant pressure). 

Eventually Tt is so high that th is smaller than the eddy 
turnover time tg = h/vc, where Vc = Vc{Tb) is the typi- 

^ The true value of j/jjo depends on the time since the last normal 
X-ray burst. This becomes important when modeling the precur- 
sor bursts. However, equating j/hc to the He ignition depth is a 
reasonable approximation when calculating the ^^C ignition condi- 
tions since the accretion timescale at yj, is ^ 10'' times longer than 
the normal burst recurrence time. 
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Fig. 1. — Evolution of the heating time = {dlnTi, / dt)~^ 
during the burst rise as a function of base temperature Tf,. Also 
shown are the thermal time at the top of the convective zone t^^ , 
the eddy turnover time at the base te, and the dynamical time at 
the base t^j. The two vertical dotted lines indicate when the burning 
makes the transition to the runaway stage and the hydrodynamic 



and an initial composition of 20% ^^C {X12 
by mass. 



0.2) and 80% ^^Fe 



cal velocity of a convective cell. At that point, convec- 
tion becomes inefficient and the fuel is rapidly consumed 
within a thin layer near the base in a local thermonuclear 
runaway. We show in § 13.21 that during the runaway, th 
becomes shorter than t^- 

Figure n shows the evolution of the heating time th — 
{dlnTb/dt)~^ during a burst rise relative to the ther- 
mal time at the top of the convective zone tth, the eddy 
turnover time at the base te, and the dynamical time 
at the base t^, for yt = 10^^ g cm~^ and X12 = 0.2. 
Here and below we define the convective, runaway, and 
hydrodynamic stages according to when t^ < th < ith, 
td < th < te, and th < td, respectively. 

3.1. Convection 

At unstable ^^C ignition, venue — 2ecooi, where 
v = dlnenuc/dlnT « 2 6 is the temperature sensi - 
tivity of the reaction fsee ICumming fc BildsteiJl2001[l . 
Thus, comparing the radiative gradient Vrad = 
(c?lnr/(iln?/)iad ~ Cnuc/ecooi to the adiabatic gradi- 
ent Vad = (c'lnr/(ilny)ad — 0.3, we have at ignition 
Vrad — 1/13 < Vad, and the burning layer is initially 
stable to convection. During the early stages of burning 
th ~ ttii and the heat flux from burning diffuses through 
the overlying accreted layers. However, the overlying lay- 
ers only heat up slightly before Vrad = Vad and a con- 
vective zone forms at the base of the burning layer. 

We account for this initial diffusive phase by comput- 
ing the one-zone rise in temperature dT^/dt = enuc/Cp 
until Vrad = Vad- To determine how the overlying lay- 
ers respond, we integrate equations ^ and |(2J) and vary 
the flux at the top boundary until the solution at the 
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Fig. 2. — Evolution of the thermal profile for an ignition col- 
umn depth of 10^^ g cm~^ and an initial composition of 20% ^^C 
and 80% ^^Fe by mass. The dashed curve shows the profile at 
ignition, the solid curve shows the profile at the onset of the con- 
vective stage, and the sequence of steep curves are the convective 
adiabats. The dot-dash curves indicate that the burning at the 
base is hydrodynamic (t^ < t^). The vertical dotted line shows 
the location of the H/He— ash interface. 



bottom boundary equals Th at yh. Figure El shows the 
change in the thermal profile between ignition (dashed 
line) and the onset of convection (solid line). Most of 
the energy release goes into heating up the immediate 
overlying layers; only a small fraction reaches the H/He 
layer (< O.li^cNo)- 

Once formed, the convective zone gradually extends 
vertically outward to lower pressures. We model the con- 
vective evo lution using a prescription similar to that of 
iWeinber'g'et al. (2006), who modeled convection during 
normal bursts. We provide a brief overview here, and 
refer the reader to their paper for further details. 

During the convective stage, we assume that the ac- 
creted layers are composed of two regions: a completely 
convective region for yc < y < yb (throughout which the 
ashes of burning are mixed) and a completely diffusive 
region for y < yc- We define the top of the convec- 
tive zone yc{Tb) as the location where the density of the 
radiative solution just exceeds that of the convective so- 
lution. We assume the thermal profile in the convective 
region follows an adiabat d\nT / d\ny = Vad(2/)- This 
is reasonable as long as the convective motions are sub- 
sonic (making convection efficient), which is true when 
tfi > te- The base temperature rises at a rate, 



dt 



Vb 

I 



■dy 



dy 



(3) 



Since the convective zone grows at the heating rate th < 
tth{yc), the radiative region cannot thermally adjust to 
the growing convective region. We therefore assume that 
heat does not diffuse out through the top of convective 
zone and the thermal profile in the radiative region is 



Fig. 3. — Minimum heating time tf^ as a function of the ig- 
nition column depth for X12 =0.1 and 0.2. The solid curves are 
models with iron as the heavy element, the dashed curves for a 
heavy composition A = 106, Z = 44. Burning becomes hydrody- 
namic for values of tf^ min that lie below the dotted curve, which 
indicates the dynamical time td at the base of the burning layer. 



unchanged from its pre-convective state. 

The thermal evolution during the convective stage is 
shown in Figure HI For yb > 10" g cm'^ and X12 = 
0.2, yc S> j/Ho when = te and the burning at the 
base undergoes a nuclear runaway well before the top 
of the convective zone ever reaches the H/He layer. It 
is therefore unlikely that convection will trigger H/He 
burning before hydrodynamic instabilities develop. 

3.2. Nuclear runaway 

As Figure^shows, by the time « 2 x 10^ K, burning 
is so rapid that th < te- Only a small fraction of the 
^^C has burned at this point and since the gas is highly 
degenerate in the burning layer, the density is nearly 
unchanged from its value at ignition. No longer able to 
transport heat via convection, the base undergoes a local 
thermonuclear runaway. 

If wc continue to naively integrate equation Q dur- 
ing the runaway, then for ignition column depths greater 
than a critical value ydym the heating time th becomes 
even shorter than the dynamical time td- In reality, a 
combustion wave must form when th ~ id- In Figure |31 
we show the dependence of the minimum heating time 
th.min on ?/{,, X12, and the composition of the burning 
layer. The nuclear runaway becomes hydrodynamical 
over the range of superburst ignition parameters. 

We can estimate how j/dyn depends on the ignition pa- 
rameters by assuming that at th = ih,min half of the 
^^C has been burned. The corresponding base tem- 
perature 16 1/2 is found by integrating dXi2/dTb « 
niir up to the midway dep letion point, with Cp 
as in lCumming- fc BildstenI 1)20011) and i?nuc the nuclear 
energy release. If we assum e th = CpTb/ enuc with enu c 
for burning to ^^Mg as in lCauehlan fc Fowled l| 1988(1 . 
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and expand the exponential in e„uc about I& 1/2, we get 
the minimum heating time in terms of the ignition pa- 
rameters, 

i,,„in « (3.5 X 10-« s) (^—) ' ' (^] . (4) 
V 5142/12/ \5X12J 

Here yu — yt/W^^ g cm^^, 314 — 5/10^'' cm s^^ and 
we assumed the equation of state is that of a rela- 
tivistic degenerate electron gas with electron fraction 
Ye. Solving th,min = td for 2/12 with td = h/cs ~ 
(7 X 10-6 s)y\!^^Y}'^ /gl{^ then gives 

0.4, (21)" 

This expression reproduces the more detailed calculation 
quite well. 

The largest thermonuclear energy release possible in a 
purely convective event E'^ax = [E^vicXi2){'iT^R^ydyn) is 

E„„., 3,4 X 10" (21)" 

(6) 

Events of less energy than this will not lead to hydrody- 
namic runaways. The observed energy release of super- 
burst s are in the range 0.5 — 1.4 x 10^^ ergs (jKuulkersI 

Hoi. 

4. SUMMARY AND DISCUSSION 

We have shown that the superburst rise evolves 
through three nuclear burning stages: an hour-long con- 
vective stage, a runaway stage, and a hydrodynamic 
stage. As rises during the convective stage, the heat- 
ing rate increases and eventually burning is so rapid that 
th ~ tf., resulting in a local nuclear runaway. During the 
runaway, th becomes even shorter than the dynamical 
time td, and a hydrodynamic combustion wave forms. 

The combustion wave propagates from the site of the 
runaway either subsonically as a deflagration or super- 
sonically as a detonation. The issues are similar to those 
found in theoretical studies of type la supernova, in 
which a C/0 white dwarf ignites -'^^C at (or near) its 
center, where the initial density is ^ 10^ g cm^'^. In 
a superburst, the combustion wave is subject to much 
larger gravitational accelerations. One can show that 



if the superburst combustion wave is a deflagration, it 
must always be subject to strong Rayleigh- Taylor insta- 
bilities. The correct solution might therefore be a deto- 
nation rather than a deflagration. 

A normal X-ray burst is found to immediately pre- 
cede the onset of the superburst in both cases where 
the RXTE Proportional Counter Array was observing 
(iStrohma ver fc Brownl l2002t iStrohmaver fc Markwarda 
2002). These precursors may be triggered by hydro- 
dynamic perturbations originating beneath the H/He 
layer. In particular, if the combustion wave propa- 
gates as a detonatio n, it will drive sho cks into the sur- 
rounding layers fsee lZingale et al.|l200lD . A deflagration 
will also generate shocks if the burning front is suffi- 
ciently distorted by Rayleigh- Taylor instabilities. The 
shocks will initially be weak, with overpressures at yh of 
^pIp ^ Enuc/Eiait ~ 0.2, where £^init is the initial in- 
ternal energy of the gas at the ignition depth. However, 
the shocks steepen as they propagate upwards, becom- 
ing quite strong by the time they reach the H/He layer 
at 1/ sa 3 X 10^ g cm~^. To estimate the strength, ap- 
proximate the initially weak, upward-propagating, shock 
as a plane parallel acoustic wave for which (Ap)'^ / pCs 
is constant by flux conservation. For y > 10 g cm"^, 
relativistic degeneracy pressure dominates and p oc y'^^^, 
Cs oc y^^^. Thus, while the shock is weak, the overpres- 
sure increases with decreasing column depth as Ap/p (x 
y-9/16^ When the shock is strong Ap/p oc y~(i+6/3)/4 
whe re - 1/(2 + \2-t/h - I 'll^^^^ (Whithamlll 971 see 
also iMatzner fc McKeelllQOl . Here 7 = 4/3 which gives 
(1 -I- 6/?)/4 = 0.56 ~ 9/16 so that over the entire range 
in y we have 

Ap/p«0.2(y/2/l2)-'/^^ (7) 

and at the H/He layer Ap/p « 20y^2^^ . In an upcoming 
paper we show that a shock of this strength deposits 
sufficient entropy into the H/He layer that it can trigger 
He burning and thereby engender a precursor to the main 
burst. 

We thank Frank Timmes for helpful conversations. 
This work was supported by the National Science Foun- 
dation under grants AST02-05956, PHY99-07949, and 
AST-0507456. 
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